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Abstract 

The distribution of the electric microfield at a charged particle moving in 
a two-component plasma is calculated. The theoretical approximations are 
obtained via the parameter integration technique and using the screened 
pair approximation for the generalized radial distribution function. It is 
shown that the two-component plasma microfield distribution shows a larger 
probability of high microfield values than the corresponding distribution of 
the commonly used OCP model. The theory is checked by quasiclassical 
molecular-dynamics simulations. For the simulations a corrected Kelbg 
pseudopotential has been used. 
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I. INTRODUCTION 



The purpose of this paper is the investigation of the microfield distribution in a two- 
component plasma at the position of a charged particle. 

The determination of the distribution of the electric microfield component created by 
one of the subsystems separately - electron or ion - is a well studied problem (for a review 
see HI]). Holtsmark f|| reduced the line shape problem to the determination of the prob- 
ability distribution of perturbing ionic electric microfield. In recent papers it was argued 
that the electric microfield low frequency part (due to the ion dynamics) also influences 
the fusion rates || and the rates for the three-body electron-ion recombination in 
dense plasmas. Holtsmark's work on the electric microfield distribution was restricted to 
ideal plasmas. The opposite limiting case of infinite coupling strength was considered by 
Mayer |5]|| within the ion sphere model. Within this model the central ion undergoes 
harmonic oscillations around the center of the negatively charged ion sphere. This results 
in a Gaussian approximation for electric microfields at the ion position. The nonideality 
of plasmas leads to quantitative corrections to Holtsmark's result as shown by Baranger 
and Mozer [0 and Iglesias |§ for the case of weakly coupled plasmas and by Iglesias et al. 
H for the case of strongly coupled plasmas. In these papers it is shown that with increas- 
ing coupling strength T the long tailed Holtsmark distribution is changed into the fast 
decaying Gaussian approximation. Here the coupling parameter T = e 2 /kTd is defined 
via the electron density n e (d = [3/47m e ] 1//3 is the average distance of the electrons). 

In the cited papers the electric microfield created by one of the subsystems has been 
studied by an almost total neglect of the influence of the other subsystem. A common 
assumption is that the distribution of the high-frequency component (due to the electron 
dynamics) is the same as that of an electron gas with uniform neutralizing background. 
This is the so called model of the one component plasma (OCP). For the ion subsystem, 
in a first approximation, the electrons are assumed to move free through the plasma. 
Since the electron motion is much more rapid than the ion one, the electrons are treated 
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as a smeared negative charged background. For simplicity this background charge was 
assumed to be uniform in the density and not to be distorted by the ion motion. This 
again is the OCP model. 

A more realistic model should also take into account the variation of the background 
charge density. A background charge distribution which differs from a uniform distribution 
results in a screening of the ion motion, the screening strength is generally frequency 
dependent, e.g. it depends on the ion velocity. In a first approximation one might neglect 
the frequency dependence of the screening. Then one arrives at the model of an OCP on 
a polarizable background (POCP). In the theory of microfields this slightly more involved 
model is used to describe the low frequency part |7|.p!0|. However, both the OCP and the 
POCP fail to describe the correlations between the electron and the ion subsystem. 

To include the electron-ion correlations one has to consider the model of a two- 
component plasma (TCP). This paper is adressed to the electric microfield studies in 
an equilibrium two-component plasma. To our knowledge the electric microfield in a 



TCP has been studied only by Yan and Ichimaru [11]. However, due to a couple of 
flaws contained in the paper of Yan and Ichimaru a further investigation is required. 
For simplicity we will restrict ourselves to the case of a two-component plasma which is 
anti-symmetric with respect to the charges (e_ = — e + ) and therefore symmetrical with 
respect to the densities (n + = rii = n_ = n e ). Further, the theoretical investigations are 
carried out for arbitrary electron ion mass ratios. To simplify the numeric investigations 
we simulated so far only a mass symmetrical (nonrelativistic) electron-positron plasma 
with m = m + = m e . We study this - so far unrealistic - case of mass - symmetrical 
plasmas in order to save computer time in particle simulations. The mass-symmetrical 
model is well suited to check the quality of various analytical approximations. In addition, 
the results of the simulation are also applicable to the case of an electron-hole plasma in 
semiconductors. 

As for the case of the OCP the microfield distribution of a TCP in the weak coupling 
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regime is approximated by the Holtsmark distribution. However, coupled plasmas are 
important objects in nature, laboratory experiments, and in technology [|T^-|T^|. 

Therefore we are interested in the modification of the microfield distribution caused 
by the coupling of plasma particles. Both theoretical investigations and semiclassical 
simulations are performed to study the microfield distribution in two-component plasmas. 

In this paper the free charges (electron and ions) are simulated by a semiclassical dy- 
namics based on effective potentials. The idea of the semiclassical method exploited in the 
numerical part of this paper is to incorporate quantum-mechanical effects (in particular 
the Heisenberg and the Pauli principle) by appropriate potentials. This method was pi- 
oneered by Kelbg, Dunn and Broyles, Deutsch and others |T5HT7| . Several investigations 
were devoted to the simulation of equilibrium two-component plasmas JT8|^22|j. Being 
interested in semiclassical methods we mention explicitely the semiclassical simulations 
of two-component plasmas performed by Norman and by Hansen JT8| , P^| . 

Certainly, such a semiclassical approach has several limits. For example, bound states 
cannot be described classically, therefore our methods are restricted to the subsystem of 
the free charges. However, this is not a very serious restriction since most of the plasma 
properties are determined by the subsystem of the free charges. 

The semiclassical approach may be very usefull to calculate a standard macroscopic 
property such as the microfield distribution since it has a well defined classical limit. The 
advantage of such an approach is the relative simplicity of the algorithm. 



II. MICROFIELD DISTRIBUTION 

Consider a two-component plasma consisting of iVj ions of one species and N e = iVj 
electrons with masses m e and mi and charges e e = — e; = — e. The total number of 
particles is N = N e + iVj. The plasma system with the total volume V and temperature 
in energy units T — 1//3 is described by the Hamilton operator 
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N a -2 1 N a N b 

H = K + V = £ + t £ £££ a6 (7v,r 6 ,) . (1) 

a=e,i i=l a a,b=e,i i=l j=l 

The interaction potential between two particles is given by the Coulomb potential 

v a b{ri,rj) = ea6b . (2) 
\ r i- r j\ 

The operator of electric field acting on a certain particle (hereafter called the first 
particle) is defined by the sum of single particle Coulomb field operators, 

N e- 
E = E E A r n) , Ej(rij) = -J-rij, th = \n -rj\ . (3) 

Define now the electric microfield distribution W(e) as the probability of measuring 
an electric field £ equal to E at the probe charge position ri, 

W(e) =< 5(e — E) > , (4) 

where < A > = (1/^) Sp (Aexp[— (3Hfj denotes the quantumstatistical average of the 
operator A, and Z = Spexp[— f3H] is the partition sum of the plasma system. 
We assume that our system is isotropic. Then we may rewrite Eq.(f|) as 

/■oo 

P(e)=2n- 1 e dl ITU) sin(d) , (5) 

JO 

where P(e) is related to W{e) by ATtW(e)e 2 de = P(e)de, and 

T(k) =< e ik - £ > (6) 

is the Fourier transform of the microfield distribution function W(e). 
It is convenient to introduce the dimensionless quantity 

F=| , (7) 

where E is defined through the total density n = N/V by E = 27r(4/15) 2 / 3 e n 2 ^ 3 . The 
probability distribution for the dimensionless field F then becomes with L = IEq, 



P(F) = 2-r- 1 F / dL LT(L) sin(FL) 
Jo 



Consider now some known limiting cases for the microfield distribution. In the weak 
coupling regime for r < 1 the Holtsmark distribution is applicable for the microfield 
distribution and we have 

T(L,r<l)=T H (L)=exp[-L 3 / 2 ] . (9) 

The other limiting case of strong coupling T 3> 1 is known so far only for the one- 
component plasma model. For the OCP the ion sphere model holds in the strong coupling 
regime. Within this model the charge will be attracted towards the center of its oppositely 
charged sphere of radius d = [3/47rn e ] 1//3 and with average density n e . The harmonic 
potential for the displacement of the center leads then to a Gaussian approximation for 
the distribution of the normalized electric field F = E/E ^ocp at the charge , 

P(F) = (2/ 7 r) 1 / 2 (F 2 /r 3 / 2 ) exp(-F 2 /2r) , (10) 

where 

4 /27r 2 \ ^ 3 

r = (bVy 1 xcothx , x = (Th 2 /4m e e 4 )T 3 , b = - I — I . (11) 

The normalizing field strength for the OCP case should be expressed in terms of the 
electron density n e only -Eo.ocp = 27r(4/15) 2 / 3 en 2 / 3 . In the case of a classical one- 
component plasma H — > the parameter r playing the role of an effective temperature in 
the Gaussian distribution Eq.fllOD simplifies and reads r c i = l/(6r), and Eq.([10|) turns 
into the expression developed by Mayer 0. 

However, there is no commonly accepted generalization of the ion sphere model for 
the two-component plasma with charges of different signs. Moreover we will show that 
in the case of TCP there is not any analogue for the Gaussian distribution in strongly 
coupled OCP. 

First we mention that the Fourier transform of the Gaussian distribution for electric 
microfield applicable in the strong coupling OCP regime equals 

LV 



T(L,r>l)=T G (L)=exp 



(12) 
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Notice that the Taylor expansion of the Gaussian function Tq{L) starts with 

Ta(L) = 1 - ^ . (13) 

On the other hand it is possible to perform exact calculations for the leading terms 
in the small L expansion of the Fourier transform T(L). From the definition of T(k) Eq. 
(H) it follows that 

T(Jfe-O) = 1 - k 2 ^^ + k 4 ^^ ±... . (14) 
y J 6 120 K J 

In Refs. fH] it was argued that it is necessary to incorporate the knowledge of the second 
moment < e 2 > into the calculation of microfield distributions in OCP. One might now try 
to generalize this idea to the case of a TCP. However, as can be easily seen the coefficient 
in the k 2 term < e 2 > diverges in the case of a TCP: 

N 

<e 2 > = < £(^(r y )) 2 +Y, E Ari3)-E h {r lh ) > . (15) 

i=2 &k 



The first sum on the r.h.s of Eq.(fL5|) can be written in terms of the partial correlation 
function of particles a and b, 

I N a N b 

<e 2 >i = g ab (r) = < Y,Y, 5 ( r ~ r ib + r ia ) > , (16) 

n a n b i = i j = i 

and reads 

4vrn e e 2 / — [g ee (r) + g ei (r)} , (17) 
Jo r z 

which diverges at small distances, since for a fluid quantum system both g e i(0) ^ and 
9ee(ty 0- m the classical OCP only g ee (r) appears with g ee (0) = 0, therefore the above 
integral and < e 2 > are finite. In contrast to what we have found Yan and Ichimaru 



|L1[ predict a finite second moment. In Ref. |TTJ no derivation of their second moment 
expression valid "strictly in the classical limit" [|ll|] is given. To isolate a possible error 
in the derivation of Yan and Ichimaru one may perform semiclassical calculations of the 
second moment. Details of the semiclassical model are given in the next Section. We 



mention here only that in the semiclassical model the quantum system is modeled by 
a system of classical particles interacting via an effective potential u a t(r) = e a eb/r + 
u s ,ab( r )> wn ere the short-range part of the effective potential u s A r ) cn ^ s the short-range 
divergency of the Coulomb potential. Therefore at short distances u s ^ a b( r — > 0) = —e a eb/r. 
Within the semiclassical model the second moment reads 

<e 2 > = ^ Gfe,(0) - g ee (0)) - /?£ h^Vu s , ab (r)) ■ 

P 0,6 X r ' 



The first term in the above equation has been reported in Ref. [fLlj| , the second term has 
been omitted. However, as may be easily seen this second term is divergent. It may be 
expressed by an integral similar to that of Eq. (|T7| ) which diverges at the lower bound. 

Thus we have established a qualitative difference between the classical OCP and the 
TCP system. For the first the second moment of the microfield distribution is finite and 
corresponds to the variance of the Gaussian distribution. 

In contrast to the OCP case the second moment of the TCP system diverges. As a 
result the TCP microfield distribution does never converge to a Gaussian distribution. 

We now generalize a coupling parameter technique which was used to calculate the 
microfield distribution of a classical OCP []8]|| to the case of a quantum TCP. Consider 
the function 

T{1) = ^ , (18) 
Z(l) = Spe ilE e~^ . (19) 

Introduce the "coupling strength" parameter A, < A < 1, of the function, 

Z(X) = Spe ixlE e- pfl . 

From the definition of T(l) in Eq. ( |T%| ) and assuming the first particle to be an electron 
one obtains 

d\nZ(X) 



lnT(/) = f 1 dX 
J o 



dX 



^2n a e a dX J dr<p(r)g ea (r,X) , (20) 



where 



il ■ r 



Z{A)Vn a n b i=lj=l 



(21) 
(22) 



The functions g a b{ r , A) may be considered as generalized partial distribution functions. 
In the case of a TCP Eq. (p0|) reads 



lnT(Z) = |e jT 1 dA | rfr0(r) [# ei (r,A) - £ ee (r,A)] 



(23) 



The above expression is still exact. The use of the "exponential approximation" (EXP) 
H ansatz leads to the expression 



g ea (r, A) ~ g ea (r, 0) exp [£ (r; A)] , a = e,i 



(24) 



with the "renormalized potential" given as 



£„(r;A) = i\l ■ E* a (r) , 

E* a (r) = E a (r la ) + dr b E b (r lb ) [g, 



Tnh) - 1 



abV'ab 



(25) 



After substitution of Eq. fl2"5D into Eq. (El) and performing the integration over A and the 
angles one gets 



T(l) 
E*Jr) 



exp 



EJr) 



An^n a / drr 2 g ea (r) ^ bo(^(r)) - 1] 
1 + 4vr^n b / r /2 c/r [# a6 (r') - 1 



E a (r) = 3 , 



(26) 

(27) 
(28) 



with jo being the Bessel function of order zero. We notice that the use of the screened 
Coulomb potential Eq.(|27|) ensures the divergency of the second moment of the TCP 
microfield distribution. In this point our theory differs essentially from the results obtained 
by Yan and Ichimaru who used a potential of mean force instead of the screened Coulomb 
field [PHI - Eqs.(p6|)-(|2"8|) constitute the so called exponential approximation (EXP) ||. It 



is known that in contrast to the so called adjustable parameter exponential approximation 
(APEX) the EXP expression poorly agrees with MD OCP data. In the APEX one 
substitutes Eq. ( PT| ) by an ad hoc ansatz for E*(r). According to this ansatz the potential 
E*(r) is approximated by a parametrized Debye potential where the parameter is choosen 
to satisfy the second moment. In order to get a generalized APEX expression for the 
TCP one should know the second moment of the TCP microfield distribution. However, 
in the above consideration we have shown that the second moment of the TCP microfield 
distribution diverges. Therefore there is not any straightforward generalization of APEX 
to the TCP case. In the weak coupling limit both approximations, EXP and APEX, 
reduce to the Debye-Hiickel (DH) approximation. 

Consider therefore the DH limit in the case of TCP. In the weak coupling limit the 
pair correlation function is given by the screened pair approximation which in our case of 
a two-component plasma reads [|T2]]: 

Pe a e b 



9a 



( r ) = S ab( r ) eX P 



(e" Kr - l) 



(29) 



where k = (47r/?X]a n a e a) 1 ^ 2 i s the inverse Debye screening length. Further 

S%\r) = const.^'exp {-(3E a ) \ ^ a | 2 (30) 

a 

is the two-paricle Slater sum written in terms of the wave functions ty a and energy levels 
E a of the pair ab, respectively. The prime at the summation sign indicates that the 
contribution of the bound states (which is not be considered here) has to be omitted. The 
Slater sum will be considered in the next section. 

To calculate the effective field E^{r) in Eq.(^) it suffices to use the linear DH approx- 
imation 

g ab (r) - 1 = - exp [-«r] , (31) 

since the nearest neighbour contribution to E*(r) is already singled out in Eq.([Z7]). In 
addition, the linear DH approximation leads to a perfect screening of the impurity charge, 
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which is an important requirement for a consistent approximation. The substitution of 
Eq.(|?TD into Eq.(P?D yields the Debye screened field 



E*(r) 



[1 + kt) exp(— Kr) 



(32) 



We put now Eqs. (|32|) and fl29|) into Eq. fl26|) and obtain the DH approximation for the 
microfield distribution in a two-component plasma. This approximation may be expressed 
in terms of the dimensionless quantities introduced in Eqs. (0) and (§) and reads 



T{L) = T ee (L) T ei (L) , 
15L 3/2 r°° dx 



In T ea (L) 



4v / 2tt Jo x 2 B(x) \ 



sin B(x) 
B(x) 



1 exp 



Z a 



Tc 



/&TL/cx 



exp (3 u s 



ex 



B{x) 



X 



1 + 



ex 





Tc ' 


j exp 





/ 2¥2 1 /3 



(33) 



( 57r )l/3 ' - e 

where the electron Wigner-Seitz radius d and the coupling constant T have been defined 



in the Introduction. Further in Eq. ( |33| ) we have introduced an effective short range 
potential 

exp(-(3u Siea (r)) = Si 2 J(r) exp (^-^^J ■ 

Equation (|8|) with T(L) from Eq. ( p3| ) constitutes the Debye-Hiickel approximation for the 
microfield distribution applicable to the weakly coupled TCP. These equations generalize 
the corresponding DH approximation used to calculate the OCP microfield distribution 
H. We mention that the approximation Eqs. ( |3"3"| ) can be directly obtained from Eq. fl2"3] ) 
using the nonlinear Debye-Hiickel approximation for the generalized radial distribution 
function, 

iMV 



9ea(r; A) 



exp 



P 



1 + 



ep 



r 



exp [-f3u s , ea (r)] 



(34) 



In the next section we consider the two-particle Slater sum and introduce the semi- 
classical model employed in the numerical simulations. 
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III. SLATER SUM, SEMICLASSICAL MODEL AND MD-SIMULATIONS 



As pointed out in the Introduction the idea of the semiclassical methods is to incor- 
porate quantum-mechanical effects (in particular the Heisenberg and the Pauli principle) 
by appropriate effective potentials. 

An easy way to arrive at effective potentials describing quantum effects is the use of 



the so-called Slater sums which were studied in detail by several authors ||i2| , |23| . The 
Slater sum caracterizes the distribution of the system in coordinate space. Choosing the 
logarithm of the Slater sum 

U^(r 1 ,...,r N ) = -ThyS(r 1 ,...,r N ) , (35) 

as a potential for the classical motion of the particles, we map our quantum system 
onto a classical one. The potentials U^(ri, . . . , r N ) are often called quantum statistical 
effective potentials and they are used to calculate the correct thermodynamic functions 
of the original quantum system fT^IJz 



The Slater sum may be considered as an analogue of the classical Boltzmann factor. 
The only modification in comparison with the classical theory is the appearance of many- 
particle interactions. If the system is not to dense (i.e., in the nondegenerate limit) one 
may neglect the contributions of higher order many-particle interactions. In this case one 
writes approximately, 

U( N \r 1 ,...,r N )n'£uiAri,r j ) , (36) 
where the effective two-particle potential u ab is defined by the two-particle Slater sum Eq. 

The Slater sum for the pair of charged particles can be approximated in different ways. 
Following Kelbg JTH] one considers the Coulombic interaction as a perturbation; in the 



first order one gets the expression 



u ab (r) = ^(F(r/X ab )), (37) 
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with 



F(x) = 1 — exp ( _ + V^ x (1 — erf (x)) 



(38) 



which we will call the Kelbg potential. Here X ab = H/ \ / 2m ab T is De Broglie wave length 
of relative motion, m~ b l = m" 1 + m^ 1 , a = e,i, m e and m 8 being the electron and ion 
masses, respectively. Further in Eq. (|37|) we have neglected the exchange contributions. 
An effective potential similar to Eq. ( |3~?D was derived by Deutsch and was used in the 
simulations by Hansen and McDonald |T9 . 



The Kelbg potential is a good approximation for the two-particle Slater sum in the case 
of small parameters ^ ab = — {e a e b ) / {TX a b) if the interparticle distance r is sufficiently large. 
At small interparticle distances it deviates from the exact value of — T • ]n(S a b(r = 0)). 
In order to describe the right behavior at small distances it is better to use a corrected 
Kelbg potential defined by |24|,|25 



u a b(r) = {e a e b /r) ■ F(r/X ab ) - k B TA ab (£ ab ) exp (-(r/A a6 ) 2 ) 



(39) 



In Eq. (^) the temperature-dependent magnitude A ab {T) is adapted in such a way 
that the Slater sum S a b(r = 0) and its first derivative S' ab (r = 0) have the exact value at 



zero distance known from previous works [12,26] . The explicit expressions read [25] 

dyyexp (-y 2 ) 



A ee = Vn\£ee\ + hi 
Aei = -V^iei + In 



exp {lc\Zee\/y) ~ 1 



V^fc(C(3) + 7 C(5)& + V< 



dyyexp (-y 2 



(40) 
(41) 



4'— 7 ' -™7 l-exp(-7r£ ei /y) 
For low temperatures 0.1 < T r < 0.3 one shall use the corrected Kelbg-potential 
Eq. (|39l) to get an appropriate approximation for the Slater sum at arbitrary distances. 
In the region of higher temperatures 

T r = T/Ti = (2Th 2 /m ie e 4 ) > 0.3 (42) 

the Kelbg potential (A ab = 0) and the corrected Kelbg potential almost coincide. At still 
higher temperatures T/Tj > 1 the Kelbg potential does not differ from the corrected Kelbg 
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potential only in the case of electron-ion interaction. For the interaction of the particles 
of the same type the correction A ab includes also the exchange effects , which make 
the potential unsymmetrical (that means u e i differ from u ee ). The potential assymetry 
becomes apparent at high temperatures (T > 100000 K) only. 

In the present work we are interested in the regime of intermediate temperatures. 
Therefore the simulations are performed with the potential Eq. ([39]) which is presented in 
Fig. 1 and compared with other potentials approximating the two-particle Slater sum. 

To check the quality of the predictions from the approximation given in Sec. [II] we 
have performed a series of molecular dynamic simulations for comparison. The leap-frog 
variant of Verlet's algorithm was used to integrate numerically the equations of motions 
obtained from the effective potential Eq. (p9|) . The simulations were performed using a 
256-particle system of electrons and positrons with periodical boundary conditions. The 
temperature of the system was choosen as T = 30 000 K, the coupling has varied from 
weak coupling (r = 0.2) up to intermediate coupling (r = 2). In the investigated range of 
plasma parameters the size of the simulation box was significantly greater than the Debye 
radius. Therefore the long-range Coulomb forces are screened inside each box and no 
special procedure like Ewald summation was implemented to calculate them. Either MD 
runs with Langevin source or MC procedures were used to establish thermal equilibrium 
in the system, both methods have led to the same results. 

In Figs. 2-5 we present the results of the approximation Eqs. (§) and (^) as well as 
the Holtsmark (Eq. (^)) approximation. The short range potential in Eq. (|33| ) is given 
by the corrected Kelbg potential without the Coulomb term 

u s , ab (r) = {e a e b /r) ■ {{F{r/X ab ) - 1)} + A ab exp(-(r/A afe ) 2 ) , (43) 

with F(x) from Eq.®. 

The results of the analytical approximation are compared with MD data. It can be 
seen from the figures that the Debye-Hiickel approximation is in good agreement with the 
MD data for the case of weak coupling, however, with increasing coupling strength this 
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agreement becomes poorer. This is not surprising, since the DH approximation is valid 
only in the weak coupling regime. To get a better agreement for the case of intermediate 
coupling one has to improve the calculation of the radial distribution function. 

From the figures we also see that the MD data show a large probability of high mi- 
crofield values. The long tails in the distribution function reflect the attraction between 
oppositely charged particles. As a result the probability to find a particle of opposite 
charge at small distances from the probe charge and thus producing large microfields is 
even higher than in the ideal Holtsmark case. This situation is in striking contrast to the 
OCP case where the repulsion of particles with the same charge leads to a small prob- 
ability of high microfield values. As for the TCP the long tails are still present in the 
case of an intermediate coupling for which the OCP microfield distribution approaches 
the Gaussian distribution Eq. ( [T0| ) |9|]. In the DH approximation the long tails are less 
pronounced for the case T = 2. Here the Debye-Htickel length is smaller than the average 
distance between the particles. Thus the particle interactions become screened even at 
short distances. A result of this unphysical screening is the supression of high microfields 
within the DH approximation and for large coupling parameters. At still higher densities 
(r > 3 at T = 30 000 K) the De-Broglie wavelength becomes comparable with the inter- 
particle distance and the semi-classical approach employed in the numerical part of the 
paper fails to describe the quantum two-component plasma properly. 

IV. CONCLUSIONS 

The electric microfield distribution at a charged particle in a two-component plasma 
has been studied. Generalizing the corresponding transformation for the case of a classical 
OCP we have expressed the Fourier transform of the electric microfield distribution in 
terms of generalized partial radial distribution functions. Using a simple Debye-Hiickel 
like generalized radial distribution function (including the unscreened short range part 
stemming from the effective potential) we have obtained theoretical predictions for the 
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electric microfield distribution of the TCP. It has been shown that in contrast to the OCP 
the second moment of the TCP microfield distribution diverges. 

Semiclassical molecular-dynamics simulations of the two-component plasma using ef- 
fective potentials have been performed . The effective potential was choosen to describe 
the nondegenrate limit of the quantum system appropriately. The microfield distribution 
for different coupling constants (from T = 0.2 to T = 2.0) has been obtained. With 
increasing coupling strength the most probable value of electric microfields is shifted to 
lower fields. However, at all coupling strengths for which the simulations have been per- 
formed the microfield distribution shows long tails indicating a large probability of high 
microfields. This behavior is in contrast to the corresponding behavior in one-component 
plasmas. It reflects the divergency of the second moment of the TCP microfield distribu- 
tion. 

At weak coupling there is an overall agreement of the microfield distribution obtained 
by the analytical approximation with the MD data. Although our simple approximation 
fails to provide accurate numerical results for larger coupling constants, the formalism 
allows to generalize the results to the case of intermediate and strong coupling. 
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FIGURE CAPTIONS 

(Figure 1) Effective potentials Eq. (|3~7D (Kelbg potential) and Eq.(^) (corrected Kelbg 
potential). The Kelbg potential is drawn for three temperatures, the corrected 
Kelbg-potential is explicitely shown at T = 10 000 K for both interactions and 
at T = 100 000 K for the electron-electron interaction only; in the other cases 
the corrected Kelbg potential coincides with the Kelbg potential within the figure 
accuracy. For comparison we have included also the low-temperature limit of the 
effective potential of free charges (the "classical" potential-dashed line); the repulsive 
part of the classical potential coincides with the bare Coulomb potential. 

(Figure 2) Comparison of microfield distribution P(F) curves at T = 30 000K and 
r = 0.2 from molecular dynamics (MD) and the analytical approximation derived 
in this work (DH) Eqs. @ and (g). 

Figure 3 Same as in Fig. 2 at T = 0.8. 

Figure 4 Same as in Fig. 2 at T = 1.2. 

Figure 5 Same as in Fig. 2 at T = 2.0. 
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